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Abstract 

The time scale for cooling of the gravitationally bound gaseous intracluster medium (ICM) is not determined by 
radiative processes alone. If the ICM is in quasi-hydrostatic equilibrium in the fixed gravitational field of the dark 
matter halo then energy losses incurred by the gravitational potential energy of the gas should also be taken into 
account. This "gravitational heating" has been known for a while using explicit solutions to the equations of motion. 
Here, we re-visit this effect by applying the virial theorem to gas in quasi-hydrostatic equilibrium in an external 
gravitational field, neglecting the gravity of the gas. For a standard NFW form of halo profiles and for a finite gas 
density, the response of the gas temperature to changes in the total energy is significantly delayed. The effective 
cooling time could be prolonged by more than an order of magnitude inside the scale radius (r s ) of the halo. Gas 
lying at a distance twice the scale radius, has negative heat capacity so that the temperature increases as a result of 
energy losses. Although external heating (e.g. by AGN activity) is still required to explain the lack of cool ICM near 
the center, the analysis here may circumvent the need for heating in farther out regions where the effective cooling 
time could be prolonged to become larger than the cluster age and also explains the increase of temperature with 
radius in these regions. 
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1. Introduction 

Clusters of galaxies are the most massive viri- 
alized objects observed in the Universe. Their po- 
tential depths correspond to virial temperatures of 
1 — 10 keV(10 7 — 10 8 K) and the baryon number den- 
sity in the inner regions could be as high as 0.1cm -3 
(e.g. Vikhlinin ct al.2005; Pointecouteau, Arnaud & 
Pratt 2005). For these temperatures and densities, 
radiative losses are expected to the bring the tem- 
perature in the central regions down to ^ 10 K 
within the available time. Yet in none of the ob- 
served clusters does the temperature drop to the 
level dictated by cooling alone. The absence of sig- 
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nificant amounts of cold gas in the cores of massive 
clusters is a major puzzle posed by X-ray observa- 
tions of massive clusters (e.g. Peterson et al.2001). 
Hence, efficient heating mechanisms must operate 
at the cores of all cooling clusters. The most popu- 
lar mechanism for suppressing cooling is energy re- 
leased by an AGN in the central cluster galaxy (cf. 
Quilis et al.2001, Babul et al.2002, Kaiser & Bin- 
ney 2003, Dalla Vecchia et al.2004, Roychowdhury 
et al.2004, Voit & Donahue 2005, Nipoti & Binney 
2005, and references therein) or by multiple AGN ac- 
tivity in all galaxies in the cores of clusters (Nusser, 
Silk & Babul 2006; Eastman et. al. 2007; Nusser & 
Silk 2008). Over- pressurized ejecta from the AGN 
transform into hot bubbles that eventually reach 
pressure equilibrium with the ICM and proceed to 
rise buoyantly away from the center. These bubbles 
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could heat the ICM by means of shock waves gen- 
erated as they expand to reach the ICM pressure 
(Nusser, Silk & Babul 2006), and by drag forces 
when they become buoyant (e.g. Churazov et. al. 
2001). Mechanical activity near the center could also 
generate sound waves which are believed to eventu- 
ally dissipate their energy in the ICM (Pringle 1989, 
Ruszkowski et al.2004, Heinz & Churazov 2005, Fu- 
jita & Suzuki 2005, Sanders & Fabian 2007). To bal- 
ance cooling in a cluster of X-ray luminosity of L x ~ 
10 44 erg s _1 , a central AGN must produce ~ 10 60 erg 
over the entire life-time of the cluster. For the most 
massive clusters (potential depths corresponding to 
velocity dispersions > 500 km/s) the required heat- 
ing could be more than an order of magnitude larger 
than the observed range of AGN energy output in 
galaxy clusters, based on the pV content of X-ray 
cavities (e.g. Best et al.2007). This is not too wor- 
rying since weak shocks could certainly compen- 
sate for the missing energy needed to balance cool- 
ing. For less massive clusters the pV energy is suffi- 
cient to balance cooling (e.g. Birzan et al.2004). The 
challenge, however, is to arrange for efficient energy 
transport from the AGN over the entire cooling core, 
or out to distances of up to ~ 100 kpc. 

The temperature in the inner regions increases 
gradually as we move away from the center. At first, 
this behavior may seem reasonable since the radia- 
tive cooling becomes more efficient nearer to the 
center. But, the cooling time is significantly shorter 
than the cluster age over a significant part of the 
inner regions and the ICM had ample opportunity 
to cool to very low temperatures (e.g. Fig. 12 in 
Wise, McNamara & Murray 2004) . So why is there 
not a temperature plateau extending over the region 
where the cooling time is shorter than the cluster 
age? One explanation might be that, on account of 
the lower density, heat conduction is more significant 
as we move away from the center. However, heat con- 
duction is not universally important in these regions 
(e.g. Wise, McNamara & Murray 2004). Here show 
that the cooling time could significantly be modified 
when the potential energy of the ICM in the dark 
halo is taken into account. We will use a version of 
the virial theorem to show that the potential en- 
ergy will absorb some of the energy loss incurred by 
the system. In some cases the potential energy will 
decrease by an amount larger than the actual loss, 
forcing the system to compensate the energy differ- 
ence by increasing its thermal energy. This is the 
case of negative heat capacity. This phenomenon is 
sometimes referred to as gravitational heating and 



has been discussed previously (e.g. Fabian & Nulsen 
1977) and is evident in numerical simulations of the 
ICM. However, the description in terms of the virial 
theorem as is done here is new and offers a simple 
analysis for assessing the dependence of the effective 
cooling time on the assumed halo profile. 

2. The virial theorem 

Hereafter we will assume spherical symmetry and 
denote by r the distance from the center. Let p g (r), 
u(r), and P = (7 — l)p u be, respectively, the gas 
density, energy per unit mass, and pressure, where 7 
is the adiabatic index. The temperature, T, is related 
to u by u — k B T/(-y — l)/m, where m is the mean 
particle mass and k B is the Boltzman constant. We 
assume a gas obeying the equation 
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where g is the gravitational force field per unit mass. 
This equation is applicable in quasi-hydrostatic 
equilibrium so that the acceleration of the gas is 
negligible. Multiplying P by r and integrating 
over the volume from r — to Ro gives the virial 
theorem, 

3(P - P )V + W = (2) 

where V = (4n/2,)B$, P = P{R ) is the exter- 
nal pressure, P = An drr 2 P(r)/V is the average 
pressure inside Ro, and the gravitational term, W, is 

Ra 

W = 4tt J r 3 p e g(r)dr (3) 


A more general derivation which includes gas mo- 
tions could be found in Ostriker & McKee (1988). 
The energy of the system in the volume V is writ- 
ten as the sum of the thermal energy jPV/ (7 — 1) oc 
Nk B T (N is the total number of particles) and the 
gravitational potential energy, U, 



E = U 



where 



PV 
~ 1 



(4) 
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and the system is assumed to reside in a static exter- 
nal gravitational potential $ and neglected gravity 
of the gas. 



U = 47r p $r dr 
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From the virial theorem ^ and the energy equa- 
tion dJ) we obtain global relations between infinites- 
imal variations (denoted by the prefix 5 ) in the to- 
tal energy, E, the thermal energy E th , V, W and U. 
Keeping a constant external pressure Pq these rela- 
tions are 



SW = 3iW - 3(7 - l)SE tb 
and 

5E = 8U + &E„ , 



(6) 



(7) 



where we have used the expression E th = PVj (7 — 1 ) 
for the thermal energy These relations must hold for 
any change in the state of the system. For radiative 
losses, the energy loss in time St is 8E = n c A(T)8t 
where n o is the electron number density and A is the 
cooling rate. Even if this energy is extracted initially 
from the thermal part, E th , subsequent evolution 
of the system will establish the relations © and 
([7]). We are working under the assumption of quasi- 
hydrostatic equilibrium so that any bulk motions 
generated during this process are neglected. In any 
case, if dissipation is important then significant gas 
motions will be converted into heat, restoring the 
above relations. 

We are set now to derive a relation between SE 
and 6E th . We write SV = {SV/SW)SW in the virial 
relation (O to obtain, 
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Writing SU = (SU/SW)SW and SV = (SV/SW)SW 
in the relation (O while taking 5W from (JSJ) we get 



SE = CSE t , 



where 



C = 1 
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The quantity C gives the ratio of the heat capac- 
ity to the standard thermodynamic heat capacity 
computed without gravity. Hence we call C the rela- 
tive heat capacity (RHC). The sign of C determines 
whether the thermal energy, E tli , and hence the tem- 
perature, T on E th /N, will increase or decrease as 
a result of changing the total energy, E. If C < 
holds, then the heat capacity is negative, i.e. the 
temperature increases when we extract energy from 
the system. For Pq = 0, the condition C < implies 

SU 1 

SW > 3(7-1) ' 



For positive RHC, C > 0, the response time of the 
gas temperature to variations in its energy is pro- 
longed by a factor C. For example, the effective cool- 
ing time is Ci cool where i cool ~ k B T/(n a A) is the 
usual radiative cooling time. 

3. applications to various forms of halo 
gravitational potentials 

We begin with the calculation of the RHC, C, for 
power-law potentials of the form, $ = A/r n so that 
g = An/r n+1 , where n / and A are constants. 
The constant A is negative for n > and positive 
otherwise. In this case W = nU and SU /SW = 1/n 
and for Pq = we have 



C = l--( 7 -l). 
n 

Thus C is negative for 
< n < 3(7- 1) , 
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which gives < n < 2 for 7 = 5/3. To estimate 
the RHS, C, for a non-vanishing external pressure, 
Pq, we need the quantity SW/5V which depends on 
gas density profile, p , in the system. We work here 
with a power-law density profile of the form, p g = 
B/r a and we compute 5W/8V under variations of 
the external radius i?o assuming that the index a 
and the total mass, Mo, inside i?o remain constant. 
Since M = in J r 2 B/r a dr we get 
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Evaluating W we get 

w = 2Qz°Lamic 



a — n 



Therefore, 



SW 
~SV 

and 

SV 
SW 



n 2 (3 



AM 



3 - a - n 4irR 3 +n 



W 
nW 



(14) 



(15) 



(16) 



(17) 



Since W < 0, this quantity is negative for n < 
and, as seen in (|10p . the heat capacity is positive 
for any Pq. For n > 0, SV /SW is positive. Thus, 
an inspection of (TIT))) reveals, for SU/SW > the 
existence of external pressure could result in C < 
even if C > for P = 0. For this to happen, the 
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value of Pq has to be adjusted such that 3PqSV/SW— 
1 is small and negative. Substituting W from the 
virial theorem ([2]) into (|17p and using the later into 
(fTOf we get the RHC in terms of the pressure ratio, 
P/Pq, as follows 



C = 1 



7 



n 1 



(18) 



n(P/P -l) 



This expression reduces to (TT21 for P/Po 3> 1. 

An intriguing case is $ = A In r and 5 = — A/r 
(A > 0). By the requirement of constant mass inside 
the varying radius i?o we find that W is constant. 
Therefore, 6U/5W is either — 00 or +00 depending 
on the sign of SU. The sign of SU is sensitive to 
the assumed density form of the gas. For example, 
for p — B/r a as above we get U = const so that 
SU = 0. For a mass distribution confined to a shell of 
negligible thickness, SU is positive when the shell is 
brought closer to the center and negative otherwise. 

We now consider the implications of the relation 
([9]) for realistic distributions of dark matter in halos. 
Therefore, we adopt the parametric form proposed 
by Navarro, Frenk & White (1996) (hereafter NFW) 
for the density profile, which is motivated by N-body 
simulations and is consistent with the distribution 
of dark matter in observed clusters (Pointecouteau, 
Arnaud & Pratt 2005). The gravitational potential 
and force field for a halo following the NFW profile 
are 



a,( r ) = _^^i n (i + cs) 
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where c is the concentration paramater, s = r/R v is 
the distance from the halo center in units of the virial 
radius R v , M v is the virial mass of the halo, and 
/ = c 2 /[ln(l + c) — c/(l + c)]. The virial mass is re- 
lated to the virial radius by M v = (47r/3)200/9 crit i? v 
where p crit is the critical cosmic density. The struc- 
ture of a halo is, therefore, determined uniquely by 
the c and R v . The scale radius r s = cR v marks the 
transition from g — const near r = to oc 1/r as 
we move further out. For this (non power-law) form 
of the gravitational field, the quantity 5U/SW de- 
pends on the assumed form of the variation in the 
gas distribution which is determined by energy gains 
and losses. We find it most instructive to focus on 
effects of local density variations on C. Therefore, 
we present here Su/Sw resulting from infinitesimal 



displacements, Sr, of a fluid element of a given mass 
as a function of its position, i.e., Su/Sw — 5<f>/5(rg). 
This type of variation is relevant for cooling/heating 
processes in a shell of finite thickness lying at a dis- 
tance ~ r from the center. We show, in figure [TJ 
Su/Sw as a function of the radius in units of the scale 
radius, r s . There a singularity at r = r th ~ 2.15r s 
which corresponds to a vanishing SW. For r > r th , 
the quantity SU/SW is positive so that according to 
(fTUj) the RHC, C, is negative (for P = 0). As we 
move to the inner regions at r < r th the quantity 
SU I SW switches signs and so the RHC becomes pos- 
itive. However, the RHC, C, is significantly larger 
than unity. At r = 0, C = 3 for Pq = and 7 = 5/3 
so that the cooling time is prolonged by a factor of 3. 
This is a modest boosting in the cooling time since 
the standard radiative cooling time scale could as 
short as 0.01 of the life-time of clusters. However, 
the dilation of cooling time scales in farther out re- 
gions could be large enough so as to exceed the clus- 
ter age. The RHC, C, is also affected by PqSV/Sw 
evaluated at the outer boundary of the shell. If the 
average pressure P inside the shell is large compared 
to Po then this term could be neglected. For realis- 
tic clusters the pressure is a steep function of radius 
for a significant part of the inner cluster cores. Of 
course the external pressure at the inner boundary 
of the shell should also be considered. However, we 
assume that the inner shell radius is small compared 
to the outer radius so that the work done by the the 
pressure on the inner shell is small. 



4. Concluding Remarks 

We used a simplified model of the ICM to study 
the its gravo-thcrmodynamical properties. In quasi- 
hydrostatic equilibrium, the inclusion of the change 
in the potential energy prolongs the response of the 
gas temperature in the inner regions lying within 
r , w 2r = 2R /c. Outside this radius, the form 

th a v / ' 

of dark halo gravitational potential is such that the 
temperature is increased as a result of energy loss, 
i.e. the gas heat capacity is negative. The boosting of 
the cooling time in the inner region formally diverges 
at r th reaching a factor of 3 as at r = 0. 

Our results may circumvent the need for heat- 
ing the ICM in regions where the standard radiative 
cooling time is an order of magnitude shorter than 
the life-time of the cluster. Those regions lie at a sig- 
nificant fraction of r s where the prolonged effective 
cooling time could be larger than the cluster age. as 
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Fig. 1. The quantity 8u/Sw = [d<t>/dr]/[d(rg)/dr] corresponding to an infinitesimal displacement, <5r, of a fluid element as a 
function of position. 



is the case for example for the cluster A1068. For 
this cluster, r s <~ 400 kpc (Pointecouteau, Arnaud 
& Pratt 2005) and the ratio of standard radiative 
cooling time to the cluster age is ~ 0.1 — 1 over the 
region between ~ 70 kpc to ~ 300 kpc. Our results 
do not eliminate the need for heating of the ICM in 
central regions (r << r th ) since the effective cooling 
time there is still shorter than the cluster age. 

As mentioned before, the approach taken here 
aims at addressing a specific point related to the 
heat capacity of the ICM. The analytic methods 
used here could be followed only by invoking a sim- 
plified (perhaps oversimplified) of the ICM. The ef- 
fect of gravitational heating has appeared in a va- 
riety of forms in the literature. For example, it has 
been discussed by Fabian & Nulsen (1977) using spe- 
cific solutions to subsonic collapse of gas in a 1/r 
gravitational potential. It is also seen in the one di- 
mensional simulations by Omukai & Nishi (1998) of 
self-gravitating gas for the formation of primordial 
protostellar clouds. It is also evident in the solutions 
for self-similar cooling flows by Bertschinger (1989). 
However, as is shown here, analysis of these prob- 
lems using virial theorem offers a simple way to un- 
derstand the physical effect and its relation to the 
halo profile and external pressure. 
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